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The amplitude for the forward electroproduction of two light vector mesons can be written com- 
pletely within perturbative QCD in the Regge limit with next-to-leading accuracy, thus providing 
the first example of a physical application of the BFKL approach at the next-to-leading order. 
Recently, a numerical determination of the amplitude has been obtained in the case of equal pho- 
ton virtualities, by using a definite representation for the amplitude and a definite optimization 
method for the perturbative series. Here, we study the main systematic effects in the previous 
determination, by considering a different representation of the amplitude and different optimiza- 
tion methods of the perturbative series. Moreover, we compare our result for the differential cross 
section at the minimum \t \ with a different approach, based on collinear kernel improvement. 
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1. Introduction 



In the BFKL approach [p, both in the leading logarithmic approximation (LLA), which means 
resummation of all terms (a s \n(s)) n , and in the next-to-leading approximation (NLA), which 
means resummation of all terms a s ((X s \n(s)) n , the (imaginary part of the) amplitude for a large-.? 
hard collision process can be written as the convolution of the Green's function of two interacting 
Reggeized gluons with the impact factors of the colliding particles (see, for example, Fig. [I]). 

The Green's function is determined through the BFKL equation. The kernel of the BFKL 
equation for singlet color representation, i.e. in the Pomeron channel, is known now both in the 
forward |Q] and in the non-forward [|3|] cases. On the other side, impact factors are known with NLA 
accuracy in a few cases: colliding partons forward jet production [Q] and forward transition 
from a virtual photon 7* to a light neutral vector meson V = p ,(0,<p [Bj. The most important 
impact factor for phenomenology, the 7* — » 7* impact factor, is calling for a rather long calculation, 
which seems to be close to completion now [0, |8|]. 

The 7* — » V forward impact factor can be used together with the NLA BFKL forward Green's 
function to build, completely within perturbative QCD and with NLA accuracy, the amplitude of 
the 7*7* — > VV reaction. This amplitude provides us with an ideal theoretical laboratory for the 
investigation of several open questions in the BFKL approach and for the comparison with different 
approaches. 

In Ref. [Q] it was shown how the 7* —* V impact factors and the BFKL Green's function 
can be put together to build up the NLA forward amplitude of the 7*7* — > VV process in the MS 
scheme and a convenient series representation for this amplitude was presented. Then, in the case of 
equal photon virtualities, i.e. in the so-called "pure" BFKL regime, a numerical study was carried 
out which led to conclude that the NLA corrections are large and of opposite sign with respect 
to the leading order and that they are dominated, at the lower energies, by the NLA correction 
from impact factors. However, a smooth behaviour of the (imaginaiy part of the) amplitude with 
the energy could be nevertheless obtained, by optimizing the choice of the energy scale so in the 
BFKL approach and of the renormalization scale }JLr which appear both in subleading terms. The 
optimization method adopted there was an adaptation of the "principle of minimum sensitivity" 
(PMS) [10] to the case where two energy parameters are present. 

Here, we want to study the main systematic effects in the determination of Ref. [Q], by consid- 
ering a different representation of the amplitude and by adopting different optimization methods of 
the perturbative series. Concerning the first effect, we consider here a representation of the NLA 
amplitude where almost all the NLA corrections coming from the kernel are exponentiated. As 
for the second effect, we compare here the PMS optimization method with two other well-known 
methods of optimization of the perturbative series, namely the fast apparent convergence (FAC) 
method JTT|] and the Brodsky-Lepage-Mackenzie (BLM) method [O]. 



Finally, we compare some of our results with those of Ref. [19], where the same process has 
been considered using some version of a collinear kernel improvement. A systematic study of the 
effect of collinear kernel improvement [13, 14, [15|, 1^, 17, 18] for the amplitude in question is in 
progress. 
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Figure 1: Schematic representation of the amplitude for the J*{p)Y'(p') ~ * V(p i)V(p2) scattering. 



2. Representations of the NLA amplitude 



The process under consideration is the production of two light vector mesons (V = p°,(0,<p) 
in the collision of two virtual photons, f*(p) f*{p') — ► V(pi) V{p2)- Here, neglecting the meson 
mass my, p\ and pi are taken as Sudakov vectors satisfying p\ = p\ = and 2{p\p2) = s; the 
virtual photon momenta are instead p = ap\ — Q\/(as)p2 and p' = a'p2 — Q\/ (ol' s)p\, so that 
the photon virtualities turn to be p 2 = —Q\ and (p') 2 = —Q 2 .. We consider the kinematics when 

1 + Q\/s + ff(s~ 2 ). In this case vector mesons 



s > Q 2 2 > A 2 QCD and a = 1 + Qj/s+< 



(s- 2 ), a' 



are produced by longitudinally polarized photons in the longitudinally polarized state [|6p. Other 
helicity amplitudes are power suppressed, with a suppression factor ~ my/Qi %. We will discuss 
here the amplitude of the forward scattering, i.e. when the transverse momenta of produced V 
mesons are zero or when the variable t = (p\ —p) 2 takes its maximal value to = —Q\Q\/s + &{s~ 2 )- 
In Ref. [[J the NLA forward amplitude has been written as a spectral decomposition on the 
basis of eigenf unctions of the LLA BFKL kernel: 



Im v 



D X D 2 



(2k) 



dv I - 



a 2 (p: R )ci(v)c 2 (v) 



l + a s (n R ) 



(v) 



+ ■ 



ci(v) c 2 (v) 



-X(v) + — + i 



dv 



+ 21n(At|) 



(2.1) 



Here, a s = CC s N c /n and D\ 2 = —47ie q fv/(N c Qi t 2), where fy is the meson dimensional coupling 
constant (f p 200 MeV) and e q should be replaced by e/v2, e/(3v2) and — e/3 for the case of 
p°, CO and <p meson production, respectively. We refer to Ref. ^ for the details of the derivation 
and for the definition of the functions of v entering this expression. Two energy scales enter 



the expression (2.1), the renormalization scale pi R and the scale so, which is an artificial scale 
introduced in the BFKL approach at the time to perform the Mellin transform from the s-space to 
the complex angular momentum plane. 
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It is easy to see that the above expression can be organized as a series: 

6162 lm s (£/ series ) 1 „ , ,2 n ~ 
D^D 2 S = j2^ aM (2 " 2) 

The b n coefficients are determined by the kernel and the impact factors in LLA, while the d„ 
coefficients depend also on the NLA corrections to the kernel and to the impact factors. For their 
expression, see Ref. [Sj. 

An alternative possibility to represent the NLA amplitude is obtained by exponentiating the 
bulk of the kernel NLA corrections, 

ImsK ) , y /,x*(».te(v)+Of(W(l(v) + ^(v)[-z(v) + ^]) 

' dv — c£(mr)ci(v)c2(v) 



D X D 2 (2%f J \s 

— oo 

1+ *wiS t S)* k (;)^(^ ttt44 1 

' (2.3) 

This form of the NLA amplitude was used in [20] (see also pi]]), without account of the last two 



terms in the second line of ( |2.3| ), for the analysis of the total 7*7* cross section. We will refer in 
the following to this representation simply as "exponentiated" amplitude. 

It is easily seen that the amplitude, in any of the given representations, is independent in the 
NLA from the choice of so and of /i« 

3. Numerical results 



In Ref. [g] we presented some numerical results for the amplitude given in Eq. ( |2.2| ) for the 
61 = 62 = 6 kinematics, i.e. in the "pure" BFKL regime. We found that the d n coefficients 
are negative and increasingly large in absolute values as the perturbative order increases, making 
evident the need of an optimization of the perturbative series. We adopted the principle of minimal 
sensitivity (PMS) by requiring the minimal sensitivity of the predictions to the change of both 
the renormalization and the energy scales, Hr and sr> We considered the amplitude for Q 2 =24 GeV 2 
and tif = 5 and studied its sensitivity to variation of the parameters /i« and Yq = ln(^o/2 2 ). We 
could see that for each value of Y = ln(s/Q 2 ) there are quite large regions in jx R and Y where the 
amplitude is practically independent on }Xr and Y and we got for the amplitude a smooth behaviour 
in Y (see the curve labeled "series - PMS" in Figs. ^| and ||). The optimal values turned out to be 
Mr — 106 an d Yo — 2, quite far from the kinematical values /J.r = Q and Yo = 0. These "unnatural" 
values probably mimic large unknown NNLA corrections. 

As an estimation of the systematic effects in our determination, we want to consider here also 



the "exponentiated" representation of the amplitude, Eq. (2.3), and different optimization methods. 1 



For more details on the following, see Ref. 
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Figure 2: lm s (,etf)Q 2 / (sD\D2) as a function of Y at Q 2 =24 GeV 2 (rif — 5): (left) series representation with 
PMS and "exponentiated" representation with PMS, (right) series representation with PMS and with FAC. 



At first, we compare the series and the "exponentiated" determinations using in both case the 
PMS method. The optimal values of pL R and Y for the "exponentiated" amplitude are quite similar 
to those obtained in the case of the series representation, with only a slight decrease of the optimal 
/is. Fig. |2| (left) shows that the two determinations are in good agreement at the lower energies, 
but deviate increasingly for large values of Y. It should be stressed, however, that the applicability 
domain of the BFKL approach is determined by the condition cc s (ij.r)Y ~ 1 and, for Q 2 =24 GeV 2 
and for the typical optimal values of one gets from this condition Y ~ 5. Around this value the 
discrepancy between the two determinations is within a few percent. 

As a second check, we changed the optimization method and applied it both to the series and 
to the "exponentiated" representation. The method considered is the fast apparent convergence 
(FAC) method [fPH, whose strategy, when applied to a usual perturbative expansion, is to fix the 
renormalization scale to the value for which the highest order correction term is exactly zero. In 
our case, the application of the FAC method requires an adaptation, for two reasons: the first is 
that we have two energy parameters in the game, /i# and Yo, the second is that, if only strict NLA 
corrections are taken, the amplitude does not depend at all on these parameters. For details about 
the application of this method, we refer to [22]. Here, we merely show the results: the FAC method 
applied to the series representation (see Fig. ^ (right)) and to the exponentiated representation (see 
Fig. |H (left)) gives results in nice agreement with those from the PMS method applied to the series 
representation, over the whole energy range considered. 

Another popular optimization method is the Brodsky-Lepage-Mackenzie (BLM) one fll2|], 
which amounts to perform a finite renormalization to a physical scheme and then to choose the 
renormalization scale in order to remove the pVdependent part. We applied this method only to the 
series representation, Eq. (22). The result is compared with the PMS method in Fig. || (right) (for 
details, see Ref. fl22|]). 

The 7*7* — > pp amplitude with the inclusion of NLA BFKL effects has been studied also in 
Ref. @. In that paper, the amplitude has been built with the following ingredients: leading-order 
impact factors for the J* —> p transition, BLM scale fixing for the running of the coupling in the 
prefactor of the amplitude (the BLM scale is found using the NLA y* — > p impact factor calculated 
in Ref. [|6|]) and renormalization-group-resummed BFKL kernel, with resummation performed on 
the LLA BFKL kernel at fixed coupling [23]. In Ref. [19] the behaviour of do jolt at t = to as a 
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Figure 4: da/dt\ t=tQ [pb/GeV 2 ] as a function of y/s at Q 2 =16 GeV 2 (n f = 5) (left) and at g 2 =4 GeV 2 
(tif = 3) (right) from the series representation with the PMS optimization method (solid lines) compared 
with the determination from the approach in Ref. [|9| (dashed lines). 



function of ->/s was determined for three values of the common photon virtuality, £2=2, 3 and 4 
GeV. 



In order to make a comparison with the findings of Ref. []19(], we computed do/dt at t = to 
for Q=2 and Q=4 GeV as functions of yfs. We used / p =216 MeV, «em = 1/137 and the two-loop 
running strong coupling corresponding to the value a s (M z ) = 0.12. The results are shown in the 
linear-log plots of Fig. ^j, which shows a large disagreement. It would be interesting to understand 
to what extent this disagreement is due to the use in Ref. [|l^] of LLA impact factors instead of the 
NLA ones or to the way the collinear improvement of the kernel is performed. 



The work of D.I. was partially supported by grants RFBR-05-02-161 1, NSh-5362.2006.2. 
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